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Abstract 

In this article I trace the historical development of tumour growth laws, which in a quantitative 
fashion describe the increase in tumour mass/volume over time. These models are usually formu- 
lated in terms of differential equations that relate the growth rate of the tumour to its current state, 
and range from the simple one-parameter exponential growth model, to more advanced models that 
contain a large number of parameters. The conclusion of this brief survey is that although much 
improvement has occured over the last century, more effort and new models are required if we are 
to understand the intricacies of tumour growth. 



The growth patterns exhibited by tumours have gathered the interest of scientists since the early 
days of cancer research, and despite almost a century of inquiry there are still uncertainties regarding 
the precise growth rate, or rather growth pattern, that solid tumours exhibit. This is in part due to 
the shift in scientific focus from the macroscopic increase in tumour size, to the (sub-) microscopic 
workings of signal transduction pathways, and the impact of specific genetic alterations. However, 
a lot can still be gained from an advance in our knowledge of tumour growth, not least in the clinic, 
where established growth curves may help to determine optimal time period between mammography 



screening (Hart et al. 1998), and may be used for calculating appropriate doses in radio- and 



chemotherapy (Steel 1989). 

The increase in size of a neoplasm, caused by the accelerated rate of cell division among the 
malignant cells, is the most basic observation to be made when studying cancer. A natural question 
to ask is then how the size of the tumour increases over time as the disease progresses. What seems 
like a very simple question, has however turned out to be immensely difficult to answer, and the 
reasons why have only become evident in the last couple of decades, when many of the complex 
processes underlying tumour progression have been discovered. 
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It was early on established that if cancer cells divide in a completely unconstrained fashion, 
and hence every cell in the tumour continuously passes through the cell cycle, and gives rise to two 
daughter cells at regular intervals, then the number of cancer cells and therefore the volume and 
mass of the tumour would increase exponentially with time. The geometric increase then implies 
that the time one has to wait for the tumour to double in size is constant over time. This picture of 
exponential growth has been shown to match the early stages of tumour growth, but in all known 
cases the doubling time eventually increases and continues to do so for the remainder of the disease. 
This can occur if either the average cell cycle time increases or if there is a loss of dividing cells 
due to quiescence or cell death. We shall discuss both these possibilities later on, but for now we 
simply note that other mechanisms need to be invoked. 

The paradigm of exponential growth is thus unable to explain growth dynamics of tumours in 
the longer term, and one has to look beyond the simple idea of unconstrained cell division to explain 
the observed data. This was indeed done in the first half of the 20th century through pioneering 



work by Mayenord (1932) and Schreck (1936) among others, who considered models where the 



growth rate is retarded. However, in order to fully appreciate this advance it is necessary to resort 
to mathematical notation. Exponential growth is described by the differential equation 



dV(t) 
dt 



rV{t) 



(1) 



which equates the rate of increase in volume dV/dt to the current volume V{t) times the growth 
rate r, assumed to be constant. This means that within an infinitesimal time interval dt the increase 
in volume dV is proportional to the current size of the tumour. The solution of this equation is 
V(t) = Voe rt , where V$ is the volume of the tumour at time t = when measurements started, and 
e is the base of the natural logarithm. Equation (JTj) can however be viewed as a special case of a 
more general equation 



dV{t) 
dt 



rV{t) 1 



(2) 



where now the rate of increase is proportional to the volume raised to the power 6, which can take 



on any value (Mendelsohn, 1963). Now if b = 1 the solution is the above exponential growth curve, 



but for b 7^ 1 the solution is given by 



V(t) = ((l-b)(rt + C)y 



(3) 



where C is a constant related to the initial condition (see fig. [T]^). By fitting growth curves from 
mouse mammary tumours to this equation, Dethlefsen et al.] ( |1968[ ) were able to show that many 
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tumours grow according to ^ with b ~ 2/3. This value of b was in fact already suggested by 



Mayenord (1932), and can be derived from simple physical considerations: Assuming that the 



tumour is spherical in shape with volume V, its surface area scales as V 2 / 3 . If we further assume 
that the growth of the tumour is limited by nutrients and/or oxygen which enter through the 
surface, then the growth rate of the tumour should be proportional to its surface area, i.e. 
In this case (b = 2/3) the solution is given by V(t) ~ i 3 , or in terms of the of the tumour radius 
R{t) = l/W) ~ t. 

In other words the radius of the tumour grows linearly with time, and this behaviour can also 
be explained by assuming that only a thin layer of cells at the surface of the tumour are in fact 



dividing. This suggestion was, at the time, highly disputed (Patt and Blackford, 1954; Baserga 



et al. 


1960; 


Mendelsohn 


1962) 



was mitotically active. In fact, for sub-exponential growth to occur, it is not sufficient for a constant 
fraction of the tumour cells to be quiescent, but an ever increasing fraction of cells must become 
mitotically inactive as the growth progresses. This seemed an unlikely scenario, given that tumours 
actually grow at a considerable rate, but we now know that this is the general mode of growth of 
solid tumours, at least prior to vascularisation: tumours contain a proliferating region of roughly 
constant width, and hence an ever diminishing fraction of active cells. 

Despite the improvement in describing tumour growth curves the above model could not account 
for the longer term reduction in growth rate often associated with later stages of the disease. 



This problem was addressed in a seminal paper by Laird (1963), who suggested a fundamentally 
different model for explaining tumour growth, namely the Gompertz model, whose unexpected 



origin deserves a mention. It was put forward by Benjamin Gompertz (1825) as a means to explain 
human mortality curves. Gompertz had made the observation that "if the average exhaustions of 
a man's power to avoid death were such that at the end of equal infinitely small intervals of time, 
he lost equal portions of his remaining power to oppose" then the number of living humans at age 
x would be given by 



L(x) 



ke 



(4) 



where k, a and b are constants, of which k and b are necessarily positive. 

The main motivation for the model was actuarial, as a practical means of determining the value 
of life insurances, and only later was it proposed as a model for biological growth. This was done 



by the geneticist Sewall Wright (1926), who observed that "the average growth power, as measured 
by the percentage rate of increase, tends to fall at a more or less uniform percentage rate" . Or in 
other words, the growth rate of an organism or organ tends to decrease at a constant rate. The 



model was first applied by Davidson (1928), for describing the growth of cattle, who also gave the 
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following derivation of the equation, which formalises what Wright had put in words. 
Assume that a growth process is governed by 



dV(t) 
dt 



r(t)V(t) 



(5) 



where V(t) is the mass or volume of an organism. This equation is similar to the equation of 
exponential growth ([!]), but with the growth rate r = r(t) now being time-dependent. We now 
assume that r(t) decreases in proportion to its current value at a constant rate p, i.e. 



dr(t) 
dt 



-pr(t). 



(6) 



The solution of these two coupled equations now yields the Gompertz equation in its more familiar 
shape 



V(t) = V Q el> 



a (1 _ e -pt) 



(7) 



where ro is the growth rate at time t = 0, and Vo is the initial tumour volume. This equation has 
successfully been fit to biological growth in a wide variety of contexts ranging the growth of internal 



organs (Laird, 1965), whole organisms (Laird et al. 1965) and entire populations (Spickett and van 



Ark, 1990). In contrast to the exponential and Mendehlson- model, the growth curve generated by 



the Gomperz equation ([7]) is sigmoidal in shape, and reaches a constant value, or asymptote, as 
t -> oo, with Voo = V e r °/ p (see fig. []]a.). 

The previous success in describing biological growth is probably what motivated the application 
of the Gompertz equation to tumour growth, but Laird also gives some justification in terms of 
tumour biology. In arguing in favour of the Gompertz model Laird notes that if only a constant 
fraction of the tumour cells were cycling then this would still give rise to exponential growth, albeit 
with a smaller growth rate. The observed exponential decrease in overall growth rate coupled with 
experimental observations suggesting that almost all cancer cells in a tumour are passing through 



the cell cycle (Laird, 1963), was instead readily explained by a model in which all cells cycled, but 



with an ever diminishing speed. 

A suggested mechanism of action was growth retarding factors, which produced by the cancer 
cells themselves, would accumulate in the growing tumour. We now know that quite the opposite is 
true; large parts of solid tumours are quiescent, and the cells that are actually dividing do this at a 



rate comparative to the one at early stages of progression (Knofel et al. , 1987). This means that the 



decaying growth rate in ^ cannot be given any natural biological meaning, since it represents a 
process which to our knowledge does not occur in real tumour. A further issue with the Gompertz 
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model is that the asympotic tumour size predicted by the model is in some cases considerably 



larger than the host animal (Laird, 1963) and also that the doubling times required to match the 



early phases of tumor growth sometimes take on unrealistically small values (< 10 hours) (Steel 



and Lamerton, 1965). suggesting fundamental problems with the modelling approach. 



The Gompertz model is however not the only model that can capture a decrease in tumour 
growth rate over time and an asymptotic mass, and at least two other models have been put 
forward as plausible candidates. The logistic (or Pearl-Verhulst) equation given by 



(It 



K 



(8) 



was first formulated by Pierre Frangois Verhulst (1838) as a means of describing the dynamics of 



a population with an intrinsic growth rate r, whose total size is limited by a carrying capacity K, 
imposed by some environmental limitation such as nutrients (see fig.[l]A.). It has since then become 
a mainstay of biomathematics and has successfully been applied to a large number of biological 



phenomena, ranging from bacterial populations to algae and mammals (Murray, 1989). Whereas 



the Gompertz equation assumes an exponentially decreasing growth rate, the logistic equation 
instead assumes that the growth rate falls off linearly with the size, until it becomes equal to zero 
when it reaches the carrying capacity V(t) = K. In terms of matching actual growth curves, the 



logistic and Gompertz equation are quite similar (Winsor 1932), with the main difference being 



that the Gompertz curve is asymmetric, with the point of inflection (the time point where the 
growth rate is maximal) occurring after 37 % of the final size has been reached, while for the 
logistic this occurs after half of the growth has occurred. The logistic equation can be derived by 
considering a spatially extended population, where reproduction is constrained by free space. But 
such a picture does not carry over to tumour growth, whose growth deceleration not only depends 
on the lack of space, but on large number of confounding factors. This leaves the model hanging 
in a phenomenological void, also inhabited by the Gompertz model, with the possible advantage 
that the parameters of the logistic have a biological motivation, describing the growth rate and the 
maximal size. 

A second candidate, which has received considerably less attention is the Bertalanffy equation 



dV ^l = aV(tflZ-bV{t) 
at 



(9) 



whose instigator was the founder of general systems theory Ludwig von Bertalanffy. It was put 



forward as a model for organism growth (Bertalanffy 1949), and its derivation is similar to the 



above model by Mendelsohn (1963), equation Q, i.e. growth occurs proportional to surface area, 
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with the additional assumption that the loss of tumour mass due to cell death occurs in proportion 
to the volume of the tumour with a constant 6, related to the commonly employed cell loss factor. 
The solution of ^ is also sigmoidal in shape, and tends to a fixed volume as time increases, 
where the growth and loss term balance each other out (see fig. [l]A.). The striking thing about 
the Bertalanffy equation is that it both matches experimental tumour growth curves well, and in 
addition has a derivation with biologically meaningful parameters. In a review of different tumour 
growth models it was in fact shown that the Bertalaffny model gave a better fit than both the 



Gompertz and Logistic model in seven out of ten cases (Vaidya and Alexandra, 1982). 

Despite this fact the Gompertz model has remained the most applied models when it comes to 



describing tumour growth curves (see for example Norton (1988)). How can this be the case? I 
believe the reasons are two-fold, one theoretical and one practical. 

The former reason is related to the nature and purpose of modeling. Do we formulate and apply 
a model in an effort to understand a system, or in order to predict its future behaviour? The answer 
is rarely clear cut, but if one leans towards prediction then a model which is disconnected from 
reality in terms of mechanisms and dynamics is acceptable, as long as it does the job of predicting. 
If one, on the other hand, one has a yearning for understanding the system at hand, then the model 
has to be derived from, and based on, real mechanisms and entities within the system. Since the 
purpose of growth curves is often to predict the future size of the tumour, it is not surprising that 
the phenomenological Gompertz model has dominated. 

The second reason is of a more practical nature, and connected to the method with which one 
fits experimental data to the Gompertz model. In order to find the parameters of the model one 
forms the quantity W(t) = log V(t)— log V(t—1) which decays exponentially with time according to 
W(t) = Ae~ pt , where A = tq/ p{e p — 1). Thus in order to find the parameters of the model one plots 
the quantity log W(t) = log A — pt as a function of time, and by applying some regression method, 
such as least squares, finds p as the (negative) slope and log A as the intercept. In this process we 
essentially take the logarithm of the tumour volume V(t) twice, suppressing any deviations from the 
Gompertz model, and thus making it quite easy to fit the model to data that follows a completely 
different growth law. This fact is illustrated in fig. [1)3, which shows four different growth curves 
(Gompertz, Logistic, Mendelsohn and Bertalaffny) and the quantities logW(t) plotted over time. 
It is evident that even though the growth curves V(t) look quite different, the plots of logW(t), at 
least for longer times, are very close to being straight lines, which is consistent with the Gompertz 
model. In the case of real tumour data with its inherent experimental error, and in some cases late 
detection, this distinction is even harder to make, and in one sense the Gompertz model is poised 
to win the race. 
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In my opinion these are the two main reasons why the Gompertz model has held a dominating 
position for nearly half a century. However, the ease with which a model can be fitted to data, 
should not determine its use over other, more accurate models, and also the latter methodological 
consideration I find worth scrutinising. A model which is grounded in the actual biology (which the 
Gompertz model is not) can not only provide predictions, but also give insight into the underlying 
dynamical process of tumour growth. By comparing the parameter values obtained from different 
patients or cell lines, one might be able to draw conclusions about the basis of growth and its 
dynamics. Among the candidates I have reviewed here it seems to me that the Bertallanffy equation 
of tumour growth is the most suitable candidate. It produces growth curves which are nearly 
indistinguishable from the well-known Gompertz model, but has the advantage of being biologically 
motivated. 

However, more research is needed in this field, which for a long time has remained dormant. A 



notable exception to this is the model by Herman et al. (2011 ), which takes into account the tumour 



vascular network and its interface with the cardiovascular system of the host. This comprehensive 
model relates tumor vascularisation and growth to metabolic rate, and gives predictions for tumour 
properties. These include growth rates, metabolic rates, degree of necrosis, blood flow rates and 
vessel sizes, all of which are clinically relevant variables, and shows that models such as this one 
are helpful in our future inquiry into the dynamics of tumour growth. 
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Figure 1: Tumour growth curves. Panel A shows the growth curves for four different models of 
tumour growth. The Gompertz (red), Bertalaffny (blue) and Logistic (black) tend to a asymptotic 
value, while the Mendelsohn model (green) gives rise to unconstrained growth. Panel B shows the 
quantity log log V(t)/V(t — 1) which is used for fitting the Gompertz model to experimental data. 
However, as the figure shows, the other models approximately exhibit a linear decrease in log W(t), 
which is usually considered to be a hallmark of Gompertzian growth. 
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